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We study itinerant ferromagnetism in a t 2 g multiorbital Hubbard system in the cubic lattice, 
which consists of three planar oriented orbital bands of d xy , d yz , and d zx . Electrons in each orbital 
band can only move within a two-dimensional plane in the three-dimensional lattice parallel to 
the corresponding orbital orientation. Electrons of different orbitals interact through the on-site 
multiorbital interactions including Hund’s coupling. The strong coupling limit is considered in 
which there are no doubly occupied orbitals but multiple on-site occupations are allowed. We show 
that in the case in which there is one and only one hole for each orbital band in each layer parallel 
to the orbital orientation, the ground state is a fully spin-polarized itinerant ferromagnetic state, 
which is unique apart from the trivial spin degeneracy. When the lattice is reduced into a single 
two-dimensional layer, the d zx and d yz bands become quasi-one-dimensional while the d xy band 
remains two-dimensional. The ground state ferromagnetism also appears in the strong coupling 
limit as a generalization of the double exchange mechanism. Possible applications to the systems of 
SrRuOs and LaAlOs/SrTiOa interface are discussed. 
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I. INTRODUCTION 

Itinerant ferromagnetism (FM) is not only a repre¬ 
sentative strong-correlation problem, but also a highly 
non-perturbative one 1-18 . It is widely known as a long¬ 
standing problem of condensed matter physics, and also 
a current research focus in ultra-cold atom physics 19-26 . 
The Stoner mechanism states that polarized electron sys¬ 
tems can save the exchange interaction energy. Never¬ 
theless, because of the associated cost of kinetic energy, 
FM is not guaranteed even in the presence of very strong 
repulsions. For example, in rigorously one-dimensional 
(ID) systems, no matter how strong the repulsive interac¬ 
tions are, the ground state is always a spin singlet, which 
is known as the famous Lieb-Mattis theorem 1 . In other 
words, electrons can remain unpolarized but avoid each 
other to reduce interaction, nevertheless, their wave func¬ 
tions are strongly-correlated. Certainly the Lieb-Mattis 
theorem in ID only applies for spin-independent systems. 
Ferromagnetism in ID is still possible if the interaction 
is spin-dependent. 

Because of the strong correlation nature of itinerant 
FM, exact theorems are important to provide reference 
points. Nagaoka’s theorem is an early example, which 
applies to the infinite U Hubbard model with a single 
hole in the half-filled background 3,27,28 . The fully po¬ 
larized FM state facilitates the hole’s coherent motion, 
which minimizes the kinetic energy of the hole and is 
therefore selected as the ground state. Another class of 
FM theorems is based on the flat band structure on line 
graphs 12-14,29,30 . Because of the divergence of density of 
states in the flat band, the kinetic energy cost because of 
spin polarization is suppressed. Metallic FM states with 
a dispersive band structure have also been proved 31 ’ 32 . 

Recently, a ground state FM theorem has been 
proved in both two-dimensional (2D) square and three¬ 


dimensional (3D) cubic lattices systems with multiorbital 
structures 33 . The band structure behaves like decoupled 
orthogonal ID chains; while, different chains are coupled 
at their crossing site through multiorbital Hubbard in¬ 
teractions. In particular, spins of each chain are not con¬ 
served but coupled by Hund’s interaction. Hence, the 
ground state FM ordering is genuinely 2D or 3D. Dif¬ 
ferent from Nagaoka’s theorem, the result of multiorbital 
FM allows a stable FM phase over a large region of filling 
factors in both 2D and 3D. An important consequence of 
this theorem is that the sign structure of the many-body 
Hamiltonian matrix leads to the absence of the quan¬ 
tum Monte-Carlo (QMC) sign problem 33 . Consequently, 
QMC simulations on finite temperature thermodynamic 
properties of itinerant FM have been performed 34 , which 
yield asymptotically exact results and shed light on the 
mechanism of magnetic phase transitions in the strong- 
coupling limit. 

In this article, we generalize Nagaoka’s theorem of itin¬ 
erant FM from the single orbital system to multiorbital 
systems. We consider the 3D cubic lattice and each site 
consists of three orbitals: d xy , d yz , and d zx . Each 
orbital has a planar orientation, and the associated band 
structure is quasi-2D like. Electrons of different orbitals 
interact through the on-site multiorbital interactions in¬ 
cluding Hund’s coupling. In the limit of intra orbital 
interaction U —> oo, states with doubly occupied orbitals 
are projected out. When each plane of the cubic lat¬ 
tice has one and only one hole in the in-plane orbital 
band, this system can be be viewed as crossing layers 
of Nagaoka FM states. We prove that, in this limit, the 
ground state of this system is the fully spin-polarized itin¬ 
erant ferromagnetic state, and it is non-degenerate apart 
from the trivial spin degeneracy. Furthermore, when this 
system is reduced to a single 2D layer system of t 2 g or¬ 
bitals, the d zx and d yz orbital bands become quasi-ID 
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and coupled to the quasi-2D band of d xy through Hund’s 
coupling. The ground state FM is still valid, where the 
quasi-ID d zx and d yz bands are allowed to take general 
values of filling, while, the d xy band can possess a single 
hole or be fully filled. Although the above exact results 
require an idealized strong-coupling limit, the strong cor¬ 
relation physics that they imply sheds important light on 
the mechanism of itinerant FM in transition metal oxides. 

The rest of this paper is organized as follows: In Sec. 
II, the multiorbital Hubbard model for the t 2 g orbital in 
the 3D cubic lattice is defined. In Sec. Ill, Theorem 1 
of the ground state itinerant FM in the 3D t 2 g orbital 
system is proved. In Sec. IV, Theorem 2 of the ground 
state itinerant FM for the reduced 2D layered system is 
proved. Discussion on connections to experiment systems 
is provided in Sec. V. Conclusions are presented in Sec. 
VI. 


II. THE MODEL HAMILTONIAN: A 3D 
MULTIORBITAL HUBBARD MODEL 

In this section, we define a 3D multiorbital Hubbard 
model in the 3D cubic lattice, which will be shown to 
possess itinerant FM ground states under conditions (I) 
and (II) in Sec. III. 

We consider a t 2 g orbital system filled with spin -1/2 
electrons; i.e., each site possesses d xy , d yz and d zx or¬ 
bitals. The Wannier wavefunction of the t 2 g orbitals is 
planar-like as shown in Fig. 1. The kinetic energy of each 
orbital band exhibits a 2D structure: Say, for electrons in 
the d xy orbital, they can only move in the zy-plane with 
a hopping amplitude tu. However, their hopping ampli¬ 
tude t± along the transverse direction of the 2 axis is very 
small. Usually, the in-plane hopping tu is assisted by the 
p orbitals of oxygen anions lying at the middle point of 
the bond, which leads to large hopping amplitudes; while, 
the transverse hopping t± can only be attributed to the 
direct overlap between two d xy orbitals offset along the z 
axis. Since d orbital Wannier functions are nearly local¬ 
ized and the z axis is perpendicular to the orbital plane, 
1 1 | is negligible in realistic transition metal oxides. Simi¬ 
larly, electrons in the d yz and d zx orbitals only hop along 
the yz and zz-planes, respectively. 

Because of the different parity eigenvalues of these 
three t 2 g orbitals with respect to the xy, yz, and zx- 
planes, they do not hybridize by the nearest neighbor 
hopping. If we neglect the longer range hopping terms, 
the kinetic energy part can simply be written as 

H K = H* + + H? x , ( 1 ) 

where H xy , H yz , and H^ x are the kinetic energies of elec¬ 
trons in the xy, yz, and zx orbital bands, respectively. 
The kinetic energy for the xy orbital band is expressed 



FIG. 1. The Wannier orbital wavefunctions of ti g orbitals: 
d xy , d yz and d zx - For electrons in the d a orbitals (a = 
xy, yz, zx), they can only move along the xy-, yz-, or zx-plane, 
respectively, but not perpendicular to the orbital orientation 
plane. 


as 

H xy = Yl h ( d lyA T ) d xyA r + a 0%) 

r 

+ d l y ,a(. r ) d xyA r + aoy) + h-c .), (2) 

where ao is the lattice constant; d xy , a (r) is the annihila¬ 
tion operator in the d xy orbital on site r with the spin- 
index cr =/ or For convenience later, we choose t\\ 
positive. For the bipartite lattice such as the cubic one, 
the sign of fy can be flipped by a gauge transformation, 
which does not affect any physical observable. The trans¬ 
verse hopping t± term is neglected in Eq. (2). Similarly, 
H y z(zx) can a ^ so defined by permuting the indices of 
orbitals and hopping directions in H xy , whose expressions 
are not repeated here. 

The interaction term is the standard multiorbital Hub¬ 
bard interaction 4 ’ 5,35 ’ 36 defined on-site as 

H 1 =U y^n a; 1 -(r)n a q(r) 

r ,a 

~ J ^2 (s a {r) ■ S b (r) - in a (r)n b (r) N j 

r,a^6 ' ' 

+ V ^2 n a {r)n b (r) 

r ,a^b 

+ A H (4,t( r ) d i,t( r )db4(r)df,,t(r) + h.c.) ,(3) 

r ,a^b 

where a = xy,yz,zx is the orbital index; n a cr (r) is the 
number of electrons occupying the a orbital at site r with 
spin-index cr; n a = n a ,t + n a, 4 .; S a (r) is the spin operator 
of the a-th orbital at site r. 

Equation (3) contains all the possible terms satisfying 
the spin SU(2) symmetry and the lattice cubic symmetry. 
The U term is the usual intra orbital Hubbard interac¬ 
tion; the V term is the inter orbital Hubbard interaction; 
the J term is Hund’s coupling with J > 0; and the A 
term describes the singlet pairing hopping process among 
different orbitals. The expressions of U, V, J, and A are 
presented in Appendix A following the standard physical 
meaning of two-body Coulomb interactions. 
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III. FERROMAGNETISM IN THE 3D t 2g 
ORBITAL SYSTEM 

In this section, we consider the 3D t 2g orbital systems 
in the cubic lattice of size L x x L y x L z . We also assume 
the following two conditions, 

(I) U -A +oo, A is finite ; 

(II) For each orbital band, there is one and only one hole 
in every layer parallel to the orbital plane. E.g., there 
is one and only one hole in every xy-plane in the d xy 
orbital band, and similarly for the d yz and d zx orbital 
bands. 

Condition II can be well-defined because of the follow¬ 
ing lemma whose proof is obvious. 

Lemma 1 The Hamiltonian of Eqs. (2) and (3) con¬ 
serves particle number in each orbital band in each plane 
parallel to the orbital orientation. 

Accordingly, the Hilbert space of the system can be fac¬ 
torized as the tensor product of the Hilbert space of each 
orbital band in each layer as 

n = 0 0 0 u%, (4) 

lz~ 1 lx~~ 1 — 1 


where, {h} represents the locations of all the holes in a 
given state and {cr} represents the spin configuration of 
this state with the labels of orbitals and layers omitted,. 
Because of the spin conservation, we can decompose the 
Hilbert space into different sectors TL Sz by the value of 
the 2 component of total spin S z , denoted as TL = © H Sz . 
Nevertheless, Tl Sz cannot be further factorized as the 
product space of different orbital bands and layers. 

Next, let us prove two lemmas as the preparation of 
the FM Theorem 1. 

Lemma 2 (Non-positivity) Under the bases 
|{A},{cr}) defined above for the Hilbert space H 
with total spin S z , the off-diagonal matrix elements of 
the many-body Hamiltonian H = Hk in + H int (see Eqs. 
(2) and (3)) are non-positive. 

Proof: : The off-diagonal matrix elements are con¬ 
tributed by the hopping part and Hund’s interaction 
part. The pairing hopping term does not exist in the 
limit of U —> +oo since states with doubly occupied or¬ 
bitals have been projected out. For the hopping term, 
because of the sign convention of the many-body bases 
defined in Eq. (6) inherited from Eq. (5), it is easy to 
check that 

({h},{a}\H t \{h'},W}) = -torO. (7) 


where l z ,x,y are the indices of the xy , yz, and 22 -planes, 
respectively. Under condition I, states with doubly occu¬ 
pied orbitals are projected out, and each orbital can only 
be occupied at most by one particle. Further, condition 
II restricts one and only one hole for each orbital band 
in a layer. In each given Hilbert space TL , (, each state 
is determined by the location of the hole and the spin 
configuration at other sites. For example, in the Hilbert 
space of the d xy orbital of the l z - th layer, we can label all 
the d xy orbitals of this layer in an arbitrary order by the 
index = 1, ■ ■ • ,L x L y . Then, the states in this layer 
can be represented as 

= ( 5 ) 

H z 


where h( v labels the location of the hole; {< 7 }^ 
represents the spin configuration; Jj[ 7 means the or¬ 
dered product of the creation operators except the 
one at the location of the hole, n- 4^fe)io) = 


4y,<7 1 ( 1 )---4 1 ^ fc _ 1 ( / C - 1 )4 y ,<7 h (/i^)4 2/ ,a fe+1 (4/ + 

l)-"4 y,a L l (L x L y ) |0) with the “hat” means the op¬ 
erator below it does not appear. Then, we can define the 
bases of the product Hilbert space for our entire system 
as 


L z L x 

IW,W) =0 01ftf'iMi.) 



L=l 


0I44W+), 

u =i 


( 6 ) 


This step is the same as that in the proof of the usual Na- 
gaoka theorem for a 2D single orbital Hubbard model 27 : 
Although here are ( L x + L y + L z ) holes in our system, the 
fermion ordering does not change under hopping because 
of Lemma 1. For the xy component of Hund’s interaction 
H Jxy = -J/2E a * b (S+S^+SfS+) with St = SZ±iS*, 
it does not change the fermion ordering either, and thus, 
its matrix elements read 

m , Ml {</}) = -J/2 or 0 , , (8) 

which are also non-positive. The V term and the 2 com¬ 
ponent of Hund’s interaction only contribute to the di¬ 
agonal part of the many-body matrix. Q.E.D. 

Let us consider a general hole and spin configuration 
satisfying conditions (I) and (II). We pick up a bond (ij) 
and consider the d a orbital of site i and the d b orbital of 
site j. If they are occupied by spin a and cr', respectively, 
let us denote this bond configuration as [d a , a (i); d btCr >(j)]. 
We have the following lemma: 

Lemma 3 The spin configuration in [d atCr (i)',db,a'(j)] 
can be flipped to [d a ,o-' (*); db,cr(j)\ by applying a series of 
hopping and Hund’s interaction processes without finally 
affecting spin and hole configurations in the rest of the 
system. 

Proof: Without loss of generality, we assume the bond 
(ij) is along the 2 axis, and only discuss how to flip 
[d a ,t(i)',d bt i(j)] to [d a ,i(i);d b ^(j)}. Since a and b can 
take any of the xy, yz, and zx, there are 9 possible or¬ 
bital configurations for a bond. Nevertheless, they can 
be classified into 4 non-equivalent classes because of the 
lattice geometry as shown in Figs. 2 (a) to (d). 
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FIG. 2. Representatives orbital configurations along a bond 
(ij) with two orbitals at sites i and j initially occupied by 
spins t and 4 respectively. Case (a): [d Z x,t(i)',d zx ,i(j)], case 
(b): [d zx ^{i)-,d yzA {j)\, case (c): [d M ,t(*); dxy,i(j)\, case (d): 

d X y : i(j)\. Any two cases of (a)-(d) are nonequivalent 
under the lattice symmetry transformation. 


For later convenience, the single hole assisted spin- 
flipping in the 2D single orbital infinite U Hubbard model 
is reviewed in Appendix B, which played an important 
role in the proof of the Nagaoka FM ground state and 
will be employed repeatedly below. 

Class (a): Let us consider a = b = zx. The same 
reasoning can also apply to the case of a = b = yz. 
Since two orbitals and the bond are coplanar and there 
is one hole in this plane, we can directly use the result in 
Appendix B to exchange their spins [d zx ^{i)\ d zx ,i(j)] —> 

[d Z x,\,{i) ) dzx, m- 

Class (&): Let us consider a = zx and b = yz , i.e., 
the configuration [d zx ^{i); d yz> ±(j)]. The reasoning be¬ 
low also applies to the case of a = yz and b = zx. Let 
us use another orbital, d yz at site i. First, we assume 
that it is occupied since we can always move an electron 
from other neighboring sites and return it back after¬ 
wards. If it is occupied by spin-4, a familiar bond con¬ 
figuration [d yz ^(i)',d yZt i(j)\ appears. As already shown 
in class (a), their spins can be exchanged to give an in¬ 
termediate configuration \d yZt \,{i)\d yz ^{j)\ for class (b). 
Then, on site i, we have both d zx ^(i) and d yz ^{i), whose 
spins can be further exchanged by the Hj term to be¬ 
come d ZXi i(i) and d yz ^{i). Combining these two steps of 
spin exchange, the initial configuration [d zX ,f (i); d yz ^{j)} 
is flipped to [d ZXt i(i); d yz ^{j)] and the third d yz (i ) orbital 
remains spin-f finally. If the d yz (i) orbital is occupied by 
spin-4, we can first apply Hund’s interaction to exchange 
spins between the d zx {i) and d yz (i) orbitals, and then ap¬ 


ply the process in class (a) to further exchange the spins 
between two d yz orbitals on sites i and j. 

Class (c) contains four equivalent configurations a = 
d zx , b d xy , a d yz , b d xy: a d xy , b d zx , 
a = d xy , b = d yz . Class ( d ) only contains one configu¬ 
ration a = b = d xy . The proof for these two classes are 
similar to that of class ( b ) by combining Hund’s interac¬ 
tion and hole’s hopping. The detailed proofs are given in 
Appendix C. Q.E.D. 

Based on Lemma 3, we can have an important prop¬ 
erty of transitivity for the many-body matrix in any sub- 
Hilbert space H Sz . 

Lemma 4 (Transitivity) Consider the Hamiltonian 
matrix in the subspace H , Sz . For any two basis vectors, 
|{/i},{u}) and !{<?},{«}}, there always exists a series of 
basis vectors \{hi}, {cti}}, \{h 2 }, { 02 }}, •••■, |{^fc},{ofc}) 
connected with nonzero matrix elements of H, such that 

({<?}, {a}\H\{h[}, W><{M, Wi}\H\ (M, W) 

x ...({h k },{a k }\H\{h},{a}) ^ 0. (9) 

Proof: Firstly, we can always apply the hopping term to 
|{/i}, {<r}) to rearrange the locations of holes of each or¬ 
bital band in each layer. Then we reach an intermediate 
state |{< 7 '}, {ck:'}) in which the locations of holes are the 
same as that in |{g}, {a}). Since the two states |{<?}, {a}) 
and |{ff , },{a / }} have the same z component of the total 
spin S z , they only differ by their spin configurations with 
a permutation of spins. 

Since any permutation can be generated by exchanges, 
it suffices to show as below that in |{g'}, (a'}) two oppo¬ 
site spins in any two orbitals can be exchanged by consec¬ 
utively applying hoppings and Hund’s interactions with¬ 
out finally affecting the configuration of the rest of the 
system. 

If the two orbitals are on the same site, it is easy to 
exchange their spins by applying Hund’s interaction with 
Hj once. If they are located at different sites, we can 
always find a path of successive bonds connecting these 
two site, and passing through nonempty sites (here are 
at most a number of min (L x , L y , L z ) sites with all three 
orbitals empty). Then, we can have a sequence of oc¬ 
cupied orbitals in which every two adjacent orbitals are 
located on two nearest neighbor sites. We can exchange 
the two spins at two ends of this path as follows: Fol¬ 
lowing Lemma 3, we can flip different spins at occupied 
orbitals on two neighboring sites. Then, by successively 
applying this operation, we are able to exchange the spins 
of two ends without affecting other parts of the system. 
Q.E.D. 

Now we are ready to prove the following theorem. 

Theorem 1 (3D FM Ground State) Consider the 
Hamiltonian Hki n + Hint satisfying conditions (/) and 
{II). The physical Hilbert space is H Sz . For any values 
ofV and J > 0, the ground states are fully spin-polarized 
and are unique apart from the trivial spin degeneracy. 
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They can be expressed as 

I^G ) = E c m,MlW’M> ( 10 ) 

where all the coefficients are strictly positive and ^2' 
means the summation over states in TL Sz . 

Proof: Because of Lemma 2 of non-positivity and 
Lemma 4 of transitivity, the Hamiltonian matrix within 
H Sz satisfies the prerequisites of the Perron-Frobenius 
theorem theorem. The importance of the transitivity to 
the non-degenerate ground state is also explained in Sec. 
Ill of the supplementary material of Ref. [33] . Then it is 
straightforward to conclude that Eq. (10) is true which 
is non-degenerate within each TL Sz . 

To show that |\t ’q z ) in Eq. (10) is a fully spin-polarized 
state, we introduce a reference state in TL Sz by summing 
over all its bases with equal weights, 

l^e/)=ElW>M)- (H) 

Since |\E'^} is symmetric under exchanging spins of any 
two orbitals, it is a fully spin-polarized state with the to¬ 
tal spin S = Ntot/2 and its z component S z . Apparently, 
('J'g l^re/) ^ Since |’Fq*} is the unique ground state 
in Ti Sz , these two non-orthogonal states must share the 
same good quantum numbers of S and S z . Q.E.D. 
Because of the spin 517(2) symmetry, all ground states 
) in different TL Sz with —N tot / 2 < S z < N tot /2 are 
degenerate, and form a set of spin multiplets with the 
maximal total spin 5 = N to t/ 2. 

Remark: Theorem 1 is true for both the periodic and 
open boundary conditions. 

Based on Theorem 1, we have the following two corol¬ 
laries with their proofs presented in Appendix D. 

Corollary 1 Under condition I and a modified condition 
II: There is one and only one doubly occupied orbital for 
each orbital band in each layer; we have that the Hamil¬ 
tonian of Eqs. (2) and (3) also possesses the fully spin- 
polarized FM ground state which is unique up to the triv¬ 
ial spin degeneracy. 

Corollary 2 If there is one and only one particle in each 
orbital band in each layer, we also have that the ground 
state is fully spin-polarized and unique up to the trivial 
spin degeneracy for any values of J > 0 and V. 

IV. FERROMAGNETISM IN THE 2D t 2g 
ORBITAL LAYER 

In this section, we will consider the same multiorbital 
Hubbard Hamiltonian of Eqs. (2) and (3) but in a sin¬ 
gle layer along the xy- plane. The d xy orbital band re¬ 
mains 2D; while, the d zx - and d yz orbitals form crossed 
ID bands with dispersion perpendicular to each other. 
The FM ground state of this system will be discussed 
when both ID and 2D bands present. 


When only the two quasi-ID bands are considered, the 
FM ground state has been proved in Ref. [33] under con¬ 
dition (I) and the following two conditions: 

(III) : Open boundary condition or periodic (anti-periodic) 
boundary condition with odd (even) number of particles 
in each row or column; 

(IV) : Arbitrary filling with at least one hole and one par¬ 
ticle in each row and each column. 

To describe the part of d zx and d yz bands with general 
fillings, let us first recapture the many-body bases con¬ 
structed for the quasi-ID system in Ref. [33] and rewrite 
them in terms of d zx and d yz bands. By Lemma 1, for any 
generic filling, we can always specify a partition of parti¬ 
cle numbers into rows X = {r, = 1, • • • , L y } and columns 
y = {ci = !,-■■ ,L X } as, M x = {N ri },Afy = {N c J, 
where N ri and N Ci are the particle numbers of d zx - and 
d yz orbitals conserved in the r»-th row and the Cj-th col¬ 
umn, respectively. We can order electrons in each row 
from the left most particle to the right most one, fol¬ 
lowed by the ordering in each column from the top to 
bottom. The corresponding many-body basis can be set 
up as 

L X Ly 

\n, s)M X ,M y =n d kc, n fi kr, i°) 

i =1 3=1 

— (j t . . . /-/t /-/t /-/t . . . /j t /71 |f)\ 

V Z 1 C L X yZ,C2 yz,Cl ZX,Vl kaj ZX,V 2 ZX,ri 1^/7 

(12) 

where j denotes the index of columns and rows; IZ = 
{rfi ; rfi | all z’s and j’s} represents the coordinates of 
occupied sites; S = {cq 3 ;/3° 3 | all *’s and j’s} represents 
their the spin configurations. The operator d\ xr . (d' yzc .) 
creates a whole line of N rj (N Cj ) d zx (d yz ) electrons in the 
row rj (column Cj) ordered from left to right (from top to 
bottom), dt x rj = lUerowr, 4x( r l)4x( r 2) ' ' ' 4z( r W,. ), 
and d' yz c . can be similarly defined. 

Now, let us consider the additional quasi-2D d xy band 
with one and only one hole. The basis for this layer of 
d xy orbital \h xy ,{a}) is defined following Eq. (5) but 
without the layer index. Then, the basis for the Hilbert 
space of this 2D system H. 2 D can be constructed by the 
direct product of the basis for the ID bands and that for 
the 2D band, 

\K,SU x ,Vy®\h xy ,{a}). (13) 

Again, because of the conservation of the z component 
of the total spin, this Hilbert space can be decomposed 
as H. 2 D = © % 2 .d ■ Following the same steps in Ref. [33] 
and in Sec. Ill, it is straightforward to show that for the 
basis defined in Eq. (13), and under condition (III) for 
ID bands, the off-diagonal matrix elements of the many- 
body Hamiltonian are non-positive. 

Below, we further show the transitivity of the Hamil¬ 
tonian matrix in the sub-Hilbert space TL^d under con¬ 
dition (IV) for d zx - and d yz bands. Since the lo¬ 
cations of electrons can be easily adjusted by apply¬ 
ing hopping terms, it suffices to show the transitivity 
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between two bases only differ by spin configurations, 
|u) = \U,S)j^ Xi j^ y ® \h xy , {ct}) and |u) = | H,S')x Xt jr y ® 

| h xy ,{a'}). Again, we only need to show that for the 
state of |it), we can exchange any two different spins by 
applying hopping and Hund’s interaction terms. If these 
two electrons are both in quasi-ID bands d zx and d yz , 
this situation has been proved in Ref. [33]. If these two 
electrons are both in the d xy band, it is reduced to the 
usual case of the Nagaoka system. 

Now let us consider the case of one electron in the 
quasi-ID bands, without loss of generality, in the d zx 
orbital band with spin-"!'; and another electron in the d xy 
band with spin-4-. We denote their locations as r zx and 
r xy , respectively. Let us identify the site r c which is in the 
same row of the d zx electron and in the same column of 
the d xy electron, and consider the d yz orbital at this site. 
We assume that there is an electron of the d yz orbital 
at r c . If not, because of condition (IV), we can always 
move a d yz electron of that column to r c by hopping. 
And the configuration in this column can be restored by 
reversing the hopping afterward. If the electron of the d yz 
orbital at r c has spin-)*, it can first be moved to r xy by 
hoppings. Then, it can exchange the spins with the d xy 
electron at r xy by Hund’s interaction. After reversing the 
hopping, this d yz electron can be moved back to r c but 
with spin-).. Further, it can be moved to r zx to exchange 
the spins with the d zx electron and be moved back to r c 
again with its original spin-f recovered. The net effect is 
the exchange of spin configurations between the d xy and 
d zx electrons without affecting other configurations. The 
case of the d yz electron at r c with spin-4. can be similarly 
proved. 

So far, we have shown both the non-positivity of off- 
diagonal matrix elements and the transitivity of the 
Hamiltonian matrix in the sub-Hilbert space Uff,. Then, 
following the same reasoning in the proof of Theorem 1, 
it is straightforward to have the following theorem 

Theorem 2 (2D FM Ground State) Consider the 
case in which there is one and only one hole in the d xy 
band. Under conditions (I), ( III ), and (IV), for any 
values of V and J > 0, the ground states are fully spin- 
polarized which is unique apart from the trivial spin de¬ 
generacy. 

Next, we consider the situation in which the d xy band 
is half-filled, i.e., there is no hole. In this case, the d xy 
band by itself is not ferromagnetic. Because of the cou¬ 
pling to the quasi-ID band, we have the following theo¬ 
rem 

Corollary 3 If the d xy band is half-filled, under the 
same conditions in Theorem 2, for any values of V and 
J > 0, the ground states are fully spin-polarized which is 
unique apart from the trivial spin degeneracy. 

Proof: We first define the basis for the local moments 
for the half-filled d xy band, which can be ordered in an 
arbitrary way as 

|{°-}> = II^^(*)I 0 >> ( 14 ) 


where (a) is an arbitrary spin distribution. Then for the 
combined system, the basis is defined as 

\H,S)M x ^ y <g) |{cr}>. (15) 

Again because of spin conservation, the Hilbert space in 
this case U' 2D can be further decomposed into the direct 
sum of different sectors of S z ’s as TL^d = ■ 

Similarly to Theorem 2, the off-diagonal elements of 
the Hamiltonian matrix is non-positive. We next show 
the transitivity of the Hamiltonian matrix in each physi¬ 
cal sub-Hilbert space H^d ■ Again, we only need to show 
that for any state in H 2 f) , opposite spins of any two elec¬ 
trons can be exchanged by applying hopping and Hund’s 
interaction without affecting other parts of the system. 
The proof is very similar to that of Theorem 2, Nev¬ 
ertheless, a new situation needs to be addressed: both 
electrons are in the d xy band with spin-f and - 4 _, respec¬ 
tively. Their locations are denoted as r and r', respec¬ 
tively. Then we can choose an electron in the d zx band, 
and, without loss of generality, assume its spin-f. Then 
according to the proof of Theorem 2, we can first flip 
the pair of electrons d zx ^ and d xy ^( r'), then their spins 
become d zx ^ and d xy ^(r'). Next, we consider the pair 
of d zx 4 and d xy ^(r) and exchange their spins. The net 
result is the exchange of the spins of two d xy electrons. 

With both results of non-positivity and transitivity, it 
is also straightforward to arrive at Corollary 3 by similar 
proof of Theorem 1. Q.E.D. 


V. DISCUSSION ON EXPERIMENTS 

Although Theorems 1 and 2 are under ideal conditions 
and limits, they do have close connections to realistic 
systems of transition-metal oxides. For the multiorbital 
Hubbard Hamiltonian of Eqs. (2) and (3), they are ac¬ 
tually a good approximation of the t 2g orbital systems 
of transition metal oxides in 3D. For example, the itin¬ 
erant FM SrRuC >3 belongs to this class of materials 37-39 , 
which is a ^-active material of 4 d electrons in a cubic 
lattice. Even though, typical interaction strength in the 
4 d electron systems are intermediately strong, it already 
exhibits the FM phase with T c = 165 K. Furthermore, 
the magnetic moment of this system is observed as 1.6/xb 
per site with the electron filling in SrRuOs as four elec¬ 
trons per site. Therefore, the FM ground state stated in 
Theorem 1 would possibly persist to the intermediate in¬ 
teraction regime and with finite electron or hole density 
away from half-filling. Nevertheless, the magnetization 
would be no longer fully polarized but partially polar¬ 
ized to save the kinetic energy cost. 

Another important system is the LaAlOa/SrTiOa in¬ 
terface between two component insulators. This interface 
is experimentally found metallic and ferromagnetic with 
large magnetization 40,41 . This is a t 2g orbital active ma¬ 
terial with 3d electrons in 2D layered systems, whose d zx 
and d yz are quasi-ID orbital bands while its d xy orbital 
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forms the quasi-2D band. For 3d electrons, the interac¬ 
tion strength is stronger than that of 4 d materials. The 
RKKY, itinerant, and double-exchange mechanisms were 
proposed to explain the FM in this system 42-4 .Here, we 
have shown that the ground state itinerant FM is fully 
spin polarized and robust for general densities in the d zx 
and d yz bands under strong intra orbital interactions. 

VI. CONCLUSIONS 

In summary, we have investigated the Nagaoka type 
itinerant FM in t 2 g orbital systems in a 3D cubic lattice. 
The hole motion in each orbital band is constrained in 
the plane parallel to the orbital orientation. Effectively, 
this system behaves as crossing planes of 2D Nagaoka FM 
coupled by on-site inter orbital Hund’s coupling. Conse¬ 
quently, 3D itinerant FM ground states are developed, 
which are proved fully polarized and unique apart from 
the trivial spin multiplet degeneracy. Also, we have con¬ 
sidered the 2D layer of t 2 g orbital systems: the quasi-ID 
bands are itinerant with arbitrary generic fillings and the 
quasi-2D band can have a single hole or be half-filled. Its 
ground state is shown remaining the fully spin-polarized 
itinerant FM. The theorems established in this article 
can be helpful for further understanding the mechanism 
of FM in SrRu 03 and the transition-metal oxides inter¬ 
face. 
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Appendix A: Expressions and physical meaning of 
U, V, J and A 

The expressions of U , V, J and A in Eq. (3) are 
standard two-body Coulomb integrals under the t 2 g or¬ 
bital basis. We assume the bare Coulomb interaction as 
V(ri — r 2 ), and express the Wannier t 2 g orbital wavefunc- 
tions (f> a { r) with a = xy,yz and zx, respectively. Then 
U, V , J and A can be represented 45,46 as 

U = J dridr 2 0*(ri)^(r 2 )U(ri - r 2 )</> a (r 2 )</> a (ri), 

V = /dridr 2C (ri)^(r 2 )^(ri-r 2 ) 

x{<M I ‘ 2 )</’ a (ri) - </> a (r 2 ) 0 b(ri)J, 

J = 2 /dridr 2 ^(ri)^(r 2 )F(n-r 2 )^(r 2 ) 0b (ri), 

A = J dridr 2 ^*(ri)^>*(r 2 )V(ri - r 2 )0 b (r 2 )<^(ri). 


where a ^ b and no summation over repeated indices is 
assumed. 

Let us explain the physical meanings of U, V, J and 
A by considering a single-site problem filled with only 
two fermions. In total there are 15 states, which can be 
classified as 3 sets of spin triplets and 6 spin singlets. The 
three sets of spin triplet states can be expressed as 

<Ai°) 

^{<T<4 + <A}|°> ( A2 ) 

with a ^ b, and their energy is V. The 6 spin singlets 
can be further classified as the orbital angular momen¬ 
tum (OAM) singlet, doublet and triplet as follows. The 
splitting between the OAM doublet and triplet states is 
because of the cubic symmetry, which is a two-particle 
analogy to the single-particle version of the f 2g and e g 
level splitting. The orbital angular momentum singlet 
state is expressed as 

^ { d ly,t d ly,l + d lz,t d lz,l + d LA*,l} 1 °)’ ( A3 ) 

whose energy is U + 2A. The orbital angular momen¬ 
tum doublet states have the energy U — A, and they are 
expressed as 

yg { d yz,t d vx,l + d lx^ d lx,i ~ 2d iy,t d iy,i} 1°)’ 

^{ d lAz,i- d Ut d l,i}\°)- ( A4 ) 

The orbital angular momentum triplet states have energy 
J + V, whose wavefunctions are expressed as 

^2 \^ d v z A dzx ,i ~ d v z ,i dzx > J [\ 

^2 \_ d ‘ zx <\ dlx v<i ~ d z x ,l dx v,t\ 

^ {<4,t4*,4 - t} i°>- ( A5 ) 

Clearly, the energy difference between the inter orbital 

singlet and triplet states is J as comes from Hund’s cou¬ 
pling. 


Appendix B: Spin flipping in a single orbital 2D 
Hubbard model in the square lattice 

To keep this paper self-contained, we review an impor¬ 
tant step showing the transitivity in the single orbital 
Nagaoka system 31 . We only consider the case of the 2D 
Hubbard model in the square lattice with U = +oo with 
a single hole [3]. The Hamiltonian can be written as 

H = tY J P{c\c j +h.c)p, (Bl) 

<y> 


(Al) 



FIG. 3. The 12 configurations with three electrons (2 spin-L 
and 1 spin-4.) and a single hole in a square plaquette. The 
hole’s motion builds up transitivity among all these states. 
From Ref. 3. 


where P is the projection operator projecting out the 
doubly occupied state. 

Consider a bond (ij) with its two sites i and j occupied 
by spins a and a' in orbitals d a and db, respectively. This 
configuration is denoted as [d a ,a(i), db, a ' (j)]- As shown 
below, the spins in this configuration can be exchanged 
to be [d a , a ' ( i ), db,<j(j)} by applying a series of hoppings in 
Eq. (Bl) without affecting hole and spin configurations 
of other sites. 

Obviously, we only need to consider the case of s z ^ s' z . 
Spin flipping can be realized by the following motion of 
the single hole. We can choose a plaquette unit contain¬ 
ing the bond (ij). If here is a hole in this plaquette, with¬ 
out loss of generality, we can assume that here are two 
spin-fs and one spin-j. in the rest 3 sites of this plaquette. 
They can form in total 12 possible combinatorial config¬ 
urations. As shown in Fig. 3, they can be connected to 
each other by simply applying hole hoppings clockwise in 
the plaquette for at most three rounds. Therefore, it is 
possible to exchange the spins on the bond (ij) without 
affecting other sites. 

If this plaquette does not contain a hole, we can first 
apply the hopping process to move the hole to this pla¬ 
quette. During this process, we require that the hole 
should not pass sites i and j, which is possible because 
even when we remove all the bonds connecting i and j, 
the remaining part of the lattice is still connected. Fol¬ 
lowing the conclusion above, we can flip the spin con¬ 
figuration on the bond (ij) without affecting other sites. 
Afterwards, we can restore the rest of the spin config¬ 
uration by reversing the hole’s motion along the same 
path on which it was brought to the plaquette before. 
Finally, the spin configuration on (ij) becomes flipped, 
i.e., [d a>IT (i),d bt<T '(j)] -A [d a y (i), d b , a (j)]. Meanwhile, 
the hole returns to its original location and spin con¬ 
figurations on other sites are restored. 


Appendix C: Exchanging spins in Classes (c) an ( d ) 

In this section, we complete the proof of Lemma 3 for 
the orbital configurations of classes (c) and (d) below. 
Proof: Class (c): We consider the case of a = zx and 
b = xy, i.e., the configuration [d zx ,\(i);d xy ^(j)\. The 
reasoning below also applies to the other 3 situations of 
a = yz, b = xy, a = xy, b = zx; and a = xy, b = yz. 
Here, the spin exchange between d zx ,f(i) and d xy ,i(j) can 
be aided by the d zx (j ) orbital. Following the reasoning 
in the main text, d zx (j ) can always be assumed occupied. 
If it has spin-t, on site j, we have d zx ^(j) and d xy ,\.(j), 
whose spins can be exchanged by Hund’s interaction to 
be d ZXt i(j) and d xy ^(j). Then bond (ij) has a new spin 
configuration [d zx ,^(i);d zx ,\.(jj\, which can be flipped as 
shown in class (a). As a result, the initial configuration 
of [d 

ZX, tW? d X y,i(j)\ is flipped to \d zx ^(i) 5 dxy, tO')] with¬ 
out affecting d zx ,^(j). Similarly, if the d zx (j) orbital is 
occupied by spin-j., we can first apply the process in class 
(a) to flip the spin configuration of d zx orbitals on sites 
i and j, and then apply Hund’s interaction to flip spins 
on d zx (j) and d xy (j) orbitals. 

Class (d): We consider the case in which both 
orbitals on (ij) are transverse, i.e., the configu¬ 
ration [d xy ^(i); d X y t i(j)]. This time we check the 
d zx (i ) orbital, and first assume it is occupied. If 
its configuration is d zx ,-\-(i), then along the bond 
(ij) we have [d zx ^(i);d xy ^(jj\, which can be flipped 
to [d zx ,\,(i); d xy ^(j)] following the steps in class 
(c). Then on site i, the spin configuration is 
changed to [d zx ,i(i); d xy ^(i)], which can be flipped to 
[d zx ^(i); d XVt i(i)] by Hund’s interaction. As a result, 
the initial configuration of [d xy ,^(i);d xy ,i(j)] is flipped 
to [d xy ,i(i);d xy ^(j)] and d zx ^(i) is maintained. If d zx (i) 
is occupied by spin-j., we can first apply Hund’s interac¬ 
tion on site i and then apply the steps presented in class 
(c). Finally, if the d zx (i) orbital is empty, we can move 
this hole to a neighboring site, and perform the above 
process, and then move the hole back. 


Appendix D: The proof of Corollaries I and II 

In this part, we prove the two corollaries in Sec. III. 

(Corollary 1) Proof: We perform a particle-hole trans¬ 
formation, i.e., d a ,a -A d\ a . Under this transformation, 
the hopping Hamiltonian Eq. (2) remains the same ex¬ 
cept for the reversed sign of t\\. Nevertheless, for the bi¬ 
partite lattice, the sign of t|| can be reversed by a gauge 
transformation, which will not change the physics. The 
physical quantities transform as follows: 

‘ft'a.a t 1 Tla,(T, S a } S a . (Dl) 

It is easy to check that for the interaction part Hi nt , U, 
V, J, and A remain the same apart from a constant and 
a term proportional to electron density. In the case of 
fixing particle numbers, the difference is just a constant 
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which does not affect real physics. Under this trans¬ 
formation, the doubly occupied orbitals are mapped to 
holes. According to Theorem 1, the ground states are 
FM states with the total spin S = N to t/ 2 — L x — L y — L z 
and are unique up to spin degeneracy. 

(Corollary 2) Proof: We order the d xy electrons layer 
by layer and define 

L z 

10*,},Ms,) = n OuU)|0>, (D2) 

lz = 1 

where r xy is the in plane location of the electron in the 
l z - th layer. Similar bases can also be defined for d yz 
and d zx electrons as | {r yz },{a} yz } and \{r zx }, {a} zx ), 
respectively. The many-body bases for the entire system 
can be defined as 

|{r}, {cr}) = I {e X y }J {<r}xy) ® I { e yz }, i 0 ”}yz ) 

<8>| {e zx },{cr} zx }, (D3) 

where {r} and {a} represent the distributions of electron 
coordinates and spins in each orbital band in each layer. 

We also need to perform a gauge transformation to flip 
the sign of f|| to be negative. Then in this case, the off- 
diagonal matrix elements of hopping are negative because 
hopping does not change the ordering of electrons in the 
definition of Eq. (D3). Because each orbital bands of 


each layer only contains one electron, only J and V terms 
contribute. Again the off-diagonal matrix elements arise 
from J, which are also negative. 

Next, we show the transitivity. Since we can also move 
the positions of electrons freely, we only need to consider 
two bases with the same electron locations but different 
spin configurations, denoted as |{r}, {er}} and |{r}, {</}). 
Then, it suffices to show that for any two electrons in 
|{r},{cr}}, we can flip their spin configuration. If these 
two electrons live in different orbitals, say, d xy and d yz , 
then the planes of their motions cross and share a com¬ 
mon line parallel to the y axis. We can move these two 
electrons to any site of this line, and then apply Hund’s 
interaction to flip their spins, and then move back to their 
original locations. If these two electrons live in the same 
orbital with opposite spins, say, two d xy electrons but in 
two parallel layers. Then we can find another electron 
in d zx orbitals. We first choose the d xy electron with 
the spin opposite to that of d zx , and switch their spins. 
Then combine the new configuration of the d zx and the 
other d xy electron, and switch their spins. The net effect 
is that two d xy electron spins are flipped, and the d zx 
electron spin is restored. 

Having proved non-positivity and transitivity, we can 
follow the same steps in Theorem 1 to prove this corol¬ 
lary, which will not be repeated here. 
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